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Abstract. We consider the survival probability of a particle in the presence of 
a finite number of diffusing traps in one dimension. Since the general solution 
for this quantity is not known when the number of traps is greater than two, we 
devise a perturbation series expansion in the diffusion constant of the particle. 
We calculate the persistence exponent associated with the particle's survival 
probability to second order and find that it is characterised by the asymmetry 
in the number of traps initially positioned on each side of the particle. 



1. Introduction 



Our understanding of the nonequilibrium dynamics of many-body systems has been 
greatly advanced through the study of particle reaction models. Typically these models 
involve a number of particles of different species that undergo reactions which lead to a 
lack of conservation of particle number. Whilst the most obvious physical application 
is to the kinetics of chemical reactions |^|, these models also enjoy mappings to 
a range of phenomena. These include interfacial growth jj), domain coarsening in 
magnetic media and fluids j|, || and aggregation ||] . Reaction-diffusion systems have 
also held their own as prototypes for the development of theoretical tools, such as the 
held theoretic renormalisation group |8| and exact methods in low dimensions @] . 

In this work we consider the dynamics of the A + B — > reaction, i.e. a system 
comprising two particle species that mutually annihilate (or form an inert product) 
on contact. This problem was introduced in the seminal paper of Toussaint and 
Wilczck as a model of monopole-antimonopole annihilation in the early universe. 
Since then, there have been applications to chemical kinetics and condensed matter 
physics — see, for example, [jfl| for a review. 



In a reaction system, such as A+B — > 0, one is chiefly interested in the concentration of 
A and B particles after time t given some initial condition. The 'traditional' approach 
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to the problem is to write down differential rate equations for the density of each 
particle species, denoted pA(x,t) and ps{x,t) respectively Q. One has for pa 

d -^=D A ^p A -rp APB (1) 
at 

and the same equation for ps if one exchanges the labels A and B. In words this 
equation says that both particle species diffuse with diffusion constants Da and Db 
and particles are removed at a rate proportional to the reaction constant r and the 
joint probability to find an A-B pair at coordinate x at time t. Under the assumption 
that the particles are well-mixed at all times, pA,B{x,t) = pA,B(t) and one can neglect 
the Laplacian in the rate equations. In physics parlance, this is equivalent to making a 
mean-field approximation which becomes exact in the "reaction-limited" regime where 
the reaction rate is slow compared to the diffusion rate. Here we are interested in the 
opposite "diffusion-limited" regime. 

The long-time form of pa and ps depends on whether the initial numbers of A and B 
particles are equal or not. Consider first the case pa(0) = /9b(0). Since each reaction 
event removes both an A and a B particle, their densities remain equal at all times. 
Then, we have for pit) — pA(t) = Pb(P) 



dp(t) _ . , 



dt 



rp Z (2) 



whose solution at large times is pit) ~ l/(rt). It turns out that this behaviour 
holds only in more than four dimensions. Through rigorous bounding arguments [ ]T2| 
and renormalisation group analysis |l3| ] it is now well understood that in dimensions 
d < 4 the density decays as p(t) ~ l/t d / i . These results serve as a reminder of the 
importance of density fluctuations and correlations in low dimensions, which were 
explicitly ignored in the mean-field approach. Also we note that, in this problem, 
there exists a dynamical scaling regime in which the characteristic length-scale in the 
problem (here the interparticle separation I = p^ x ^ d ) is related to a characteristic 
time-scale through a power-law. 

When the initial densities of A and B particles are not equal, the situation is somewhat 
different. Consider an initial condition that has pa(0) < Ps(0). Then, as time 
progresses the A + B — > recation implies that the ratio p A if) I PB(t) decreases. 
Eventually one has a small number of A particles in a sea of B particles whose density 
remains effectively constant. In the rate-equation formalism (|IJ) we have, under the 
assumption that particles are well-mixed and PB{i) = const — pb, that the density of 
A particles decays as ~ exp(— rpst). Again it is known that this result holds 

only in a space of suitably large dimensionality. In particular Bramson and Lebowitz 
|fl2[ used bounding arguments to show that, at long times, the A particles experience a 
stretched-exponential decay PA(t) ~ ex p( — A<j^ 2 ) in fewer than two dimensions, and 
PA(t) ~ exp(— \2t/ lni) in the critical dimension d = 2. Values of the constants 
for d < 2 were recently reported by us juj with the details to be presented elsewhere 



Clearly the special case of a single A particle in the presence of a sea of B particles will 
be governed by these asymptotics. Then, one can interpret the p A if) as the survival 
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probability of the A particle. If the distribution of B particles is homogeneous over all 
space, and the diffusion constants of the A and B particles are the same, one can also 
view PA(t) as the fraction of particles that have not met any other particles. Thus 
the reaction A + B — > in the limit of a low density of A particles has been discussed 
under the guises of uninfected walkers ]l6| in which random walkers infect each other 
on contact, diffusion in the presence of traps |l7], [T^] in which the B particles are 
considered as traps for the A particles, and predator-prey models |H| in which one 
asks for the survival of a prey (the A particle) being 'chased' by diffusing predators 
(the B particles). To avoid confusion, we shall use only the trapping terminology, 
and in this work we restrict ourselves to the case of a single A particle and a finite 
number of diffusing traps (B particles) in one spatial dimension. When we make no 
distinction between the particle and the traps, we shall use the generic term walkers. 

Before we formally introduce our model and review some exactly known results, let 
us make a few comments about the relationship between the trapping reaction and 
persistence in diffusive systems. Persistence is a property that can be defined for any 
stochastic dynamical system as it is simply the probability that a random variable does 
not change sign between time and a time t. In many cases, a persistence probability 
Q(t) decreases as a power-law in time, Q(t) ~ t~ 9 in which 6 is a persistence exponent 
(see, e.g., |2(| for a review). In the trapping reaction we can define two distinct 
persistence properties: (i) a site persistence, which measures the probability that a 
particular point in space (or site on a lattice) has not been crossed by any walker; 
and (ii) a walker persistence, which is the probability that a particular walker (here 
the A particle) has not yet met another walker. In this terminology, one can make 
the transition from walker to site persistence by decreasing the particle's diffusion 
constant to zero. Then the particle is static, and the probability that it survives is 
equal to the probability that a particular point has not been crossed. 

In the following section we will show that the site persistence probability in the 
presence of N traps in one dimension is known exactly to decay with an exponent 
9 = N/2. On the other hand the more general walker persistence probability, where 
the particle has a nonzero diffusion constant, is not known exactly. The principal 
purpose of this paper is present the full details of a perturbation theory for the walker 
persistent exponent that was outlined in a previous work [jl4|. As well as describing 
the approach in greater detail, we also extend the calculation to second order in 
(a quantity closely related to) the particle's diffusion constant. We find that the 
persistence exponent depends on a quantity that characterises the level of asymmetry 
between the number of traps initially positioned to each side of the particle. 

To close this introductory section, we make a few remarks as to why the case where 
both the particle and the traps diffuse is considerably more complicated than that 
in which either species is static (the case of static traps is sometimes referred to 
as the Donsker-Varadhan problem jL7) for which the walker persistence decays as a 
stretched exponential) . When the walker is static, the probability that none of the N 
independently diffusing traps has crossed a particular point is simply the product of 
probabilities for each of the traps separately not to have crossed that point. However, 
any movement of the particle (in the subset of trajectories for which the particle 
survives) correlates the motion of the traps, and so the probabilities for each of the 
traps not to have met the particle are no longer independent. To understand this 
point more clearly, consider a set of walkers on the one-dimensional lattice. When the 
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particle takes a step to the left (say), all the traps to the left hop one site closer to the 
particle in concert, whereas those to the right all hop one site away. It is this collective 
motion of the traps in the particle's frame of reference that makes the problem very 
difficult mathematically, and as far as we are aware, exact expressions for the walker 
persistence exponent are known only when N < 2. 

In the following section, we discuss these known results more quantitatively, and 
make the connection to isotropic diffusion in an iV-dimensional 'hyperwedge' geometry. 
Then, in section |^, we describe how one sets up a perturbation theory using a path- 
integral approach. A judicious choice of timescale reveals that the natural parameter in 
the perturbation expansion has a geometric significance in the hyperwedge problem. In 
section |J we go on to calculate the coefficients of the first and second order terms in the 
perturbation expansion, before concluding with some discussion and open questions. 

2. Model definition and review of known results 

The model we consider is defined as follows. At time t = 0, a particle is placed at the 
origin of a one-dimensional space Xo(0) — 0. Additionally, a set of mobile traps are 
placed at positions Xi(0) — yi with i = 1, 2, . . . , iV such that N is the number of traps. 
We will at times need to distinguish between those traps initially placed to the left 
and right of the particle. To this end we introduce the shorthand Ui = sgn(j/j) and the 
quantities Nl and Nr which are the number of particles that have <x; < and crj > 
respectively. 

The system evolves through the independent diffusion of the particle and traps. We 
shall take the diffusion constant of the particle to be D' and that of the traps to be 
D. Then, we can write a Langevin equation 

Xi = rn{t) (3) 
in which T]i(t) is a Gaussian white noise with zero mean and correlator 

{r H (t)Vj(t'))=2D i S i!j 6(t'-t) , (4) 
where D = D' and D; = D for 0. 

If, at any time, a trap coordinate Xi(t) coincides with the particle coordinate xo(t), the 
particle becomes 'trapped'. We seek an expression for the survival probability Q(t), 
i.e. the probability that the particle has not yet fallen into a trap. 

It is convenient to work in the particle's frame of reference and so we introduce the 
relative coordinates Xi — Xi — xq an d Yi = yi- Then the connection to persistence 
becomes clearer, since Q(t) is the probability that none of the variables Xi has changed 
sign until a time t. From (^) and (||) we obtain the underlying Langevin equations in 
the variables Xi. These read 



(5) 



Perturbation theory for the one- dimensional trapping reaction 



5 



in which the noise has mean zero and correlator 

(60)&(*')) = 2{P&ij + D')S(t' -t) = 2Dij5(t' - t) . (6) 

Here we have introduced the matrix = DSij + D' which expresses the way in 
which the trap's trajectories are correlated in the frame of reference of a diffusing 
particle (i.e. when D' ^ 0). 

In the introduction, we remarked that the survival probability Q(t) can be calculated 
exactly when D' = 0. Although this is a classic result (see, e.g., jyj), we present the 
details of the calculation here as it is will be needed when calculating the coefficients 
in the perturbation expansion that we construct for general D' in the next section. 

We begin by writing down the Fokker-Planck equation for the probability distribution 
function G N {{X},t\{Y}) that holds when D' = 0. It reads 

^G N ({X},t\{Y})=Dj2-^ ? G N ({X},t\{Y}) (7) 

and is subject to the absorbing boundary condition Gat = if any of the X; are 
zero (i.e., when a trap meets the stationary particle). Clearly, since the traps diffuse 
independently, the solution of this equation is just 

N 

G N ({X},t\{Y}) = l[G 1 (X i ,t\Y i ) (8) 

i=l 

where Gi(X, t; Y) is the probability for a walker starting at Y to be at X and not to 
have crossed the origin up to time t. For a particle starting at Y > 0, this quantity is 
given by 



Gx(X,t\Y) 



1 



\jAitDt 



(X-Y) 2 \ ( (X + Y) 



exp — exp 



ADt 1 V 4L>i 



(9) 



for X > 0. We note that this is the correct solution since the two terms on the right- 
hand side separately obey (Q) with N = 1 and this particular combination of solutions 
satisfies the boundary condition Gi(0,t\Y) = 0. 

To obtain the probability for a walker starting at Y not to have crossed the origin 
until time t we simply integrate (||) over all X > to find 

Qi(t\Y) = erf . (10) 



Note that, by symmetry, we must have Qi(t\ —Y) = Qi(t\Y) so the general result is 
Qi(t|Y) = exf (J^) . (11) 



2VDt 
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Finally, using (g) we find the probability that none of the traps has met the particle 
is 

' ' IKI \ /in"' 2 " 



in which we have used the fact that 

« fa = ^(*-T + -) (13) 

to obtain a form for Qjsr(t) valid at large times. The persistence exponent 9 associated 
with the survival probability Q(t) can be defined formally as 

9 = - lim t^-\n.Q(t) (14) 
t^cc dt 

and so for the case D' = we have from (|l2| ) that 9 — N/2. 



It is instructive to consider an approach that has been used (see 19, |22J for a more 
in-depth overview) to find the persistence exponent 9 exactly for the case N = 2 and 
D' =/= since it will reveal a connection with the perturbation series that we derive in 
the next section. 

The essence of this approach is to perform a second coordinate transformation that 
renders the correlated diffusion described by equations (|5|) and (j^) isotropic in the 
TV-dimensional space spanned by the coordinates To perform this transformation 
one considers the diffusion matrix = DSij + D' . The eigenvectors of this matrix 
indicate the directions in which the iV-dimensional diffusion is independent whilst 
the corresponding eigenvalues give the diffusion constants in those directions. One 
finds that one of the eigenvectors is (1, 1, 1, . . . , 1) and has an eigenvalue D + ND'. 
The remaining eigenvectors are degenerate with eigenvalue D and hence form a basis 
set in the (N — l)-dimensional space perpendicular to (1, 1, 1, . . . , 1). In order to 
make the diffusion isotropic, one must scale in the (1, 1, 1, . . . , 1) direction by a factor 
(1 + ND'/D)- 1 ' 2 . 

Whilst this coordinate transformation gives rise to simple, isotropic diffusion in 
the iV-dimensional space, the boundary conditions become more complicated. In 
the space spanned by the variables Xi, the absorbing boundary is constructed 
from the set of orthogonal planes that have Xi = with i — 1,2,..., TV. The 
coordinate transformation just described rotates these planes so that they are no 
longer orthogonal. It is a straightforward exercise in linear algebra to determine 
that the angle between two planes i and j in the transformed coordinate system is 
&ij = arccos(— (TiCjaAr) in which the parameter ajv takes the form 

aN = D + {N-1)D> ■ (15) 

Hence two planes that correspond to traps that start on opposite sides of the particle 
close up whereas two planes that correspond to traps starting on the same side of the 
particle open out. 
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To determine the persistence exponent 9 given by equation (|14|) it is necessary to 
solve the diffusion equation in an iV-dimensional wedge geometry. As far as we know, 
results are available only for the case of a two-dimensional wedge |2^, ^2|. Then, the 
particle's survival probability decays with the exponent 

(16) 



26 12 

For the case where both traps initially surround the particle, we have o~\ct2 = — 1 and 
so 

(17) 



2 arccos a 2 

Similarly when both traps are initially positioned to one side of the particle, we have 
ci<7 2 = 1 and thus 

(18) 



2(7r — arccos a 2 ) 
since arccos(— x) = 7r — arccos(x). 

Although this approach has thus far proved intractable for N > 2, we have spent some 
time outlining it because, in the perturbation theory that we construct below, it turns 
out that the natural expansion parameter is the quantity ajv introduced above. Note 
that ajy = corresponds to the case where the particle is static and the A-dimensional 
hyperwedge is in fact a corner of an A-dimensional hypercube. Thus the expansion 
for small derived below using an alternative path-integral approach gives results 
for diffusion in a geometry perturbed from an A-dimensional hypercube. 



3. Perturbation expansion via a path-integral approach 



The starting point in the derivation of a perturbation expansion for the persistence 
exponent 8 is to note that the statistical weight a particular set of trajectories 
X(t) = [Xx(t),X2(t), . . . ,Xjv(i)] governed by equations (||) and @ can be written 
in the form exp(— S*[A]). In this expression, S[X] is an 'action' functional defined as 

i N r* 

S[X] = -J2 [D-% / dt' X^Xjit') (19) 

i,3'=l 

in which D~ x is the inverse of the matrix D in equation (^|). It is straightforward to 
show that 

" 5 ( 4 * - d^Wd) ■ < 20 > 

It is helpful at this stage to introduce a rescaled time variable 
D + ND' 



D + (N-1)D' 



Dt (21) 
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such that the action becomes 

D' 

S[X] =S - D + ^ N _ ^ D/ Sl = S o~ a N Si (22) 

in which the diagonal part Sq and off-diagonal part Si of the action are given by 

So[I] = jE [at'IUt')] 2 (23) 

S ^ = lJ2 fdr'X^X^r') (24) 

#j 

where now a dot denotes differentiation to the rescaled time variable r. Note that 
when a at = the action S*[X] = 5o[^] which gives the statistical weight for TV 
walkers executing independent diffusive motion with diffusion constant unity. We will 
treat the off-diagonal part (that describes the correlations between the walkers) as a 
perturbation with the quantity a at given by equation ( fl5|) as the expansion parameter, 
as advertised. 

To this end, observe that out of all possible trajectories X(t), only a subset 
corresponds to the particle surviving until time r. These trajectories are defined 
through the persistence condition that sgnX(r') = sgnX(O) for < r' < r. We shall 
use the symbol J r VX(t) to denote integration over this restricted set of surviving 
trajectories. Then, the probability that the particle survives is 

= J R VX(r)eMS[X]) 
fVX(T)exp(-S[X}) 

in which the integral of the weights over the restricted trajectories is normalised by 
the integral over all possible trajectories. 

Now, note that the survival probability Qn(t) of a stationary particle in the presence 
of N traps diffusing independently with diffusion constant unity is given by the path- 
integral expression 

= /j? ^(,)exp(-5 W) 
JVX(r)e X p(S [X}) 



Combining this with ( |25[ ) we find that 

0(t) = q (t) S R VX{r)eM~S[X]) J VX(r) exp(-S [X]) 
J R VX(T)eM-S [X}) JVX(T)eMS[X}) 

= Qn(t) ) s > R (27) 

in which the notation (•) indicates an average over paths weighted by exp(— Sq[X]) 
and the subscripts R and U denote the restricted and unrestricted ensembles of paths. 
(Recall that paths that cross the origin are excluded from the former). 
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To obtain a perturbative expression for the persistence exponent 9 we make use of 
(H) and (|27|) to find, in the rescaled time variable r, 



lim t— In Q(t) 
t^oo dr 



lim t— 

t— >oo dr 



(lnQ A ,(r)+ln(e Q « Sl ^]) fl -ln(e tt "' Sl ^]) c/ ) • (28) 

The asymptotic form of Qn{t) is known from equation (|l2|). To compute the averages 
of perturbatively, we make use of the fact that ln(exp Ax) defines the generating 

function of the cumulants of x. Specifically 

ln(e A *) = \{x) + y ((x 2 ) - (x) 2 ) + O (A 3 ) . (29) 



Putting all this together, we find that 8 has a series expansion 

(30) 



N d 
— — lim r — 

2 r^oo dr 



a w ((Si)* - {Sx) v ) + ((Sl) R - (Si)% - + 



to second order in the parameter defined in equation ()l5|). Of course, one could 
go to higher order by including more terms from the cumulant expansion (E9). 



4. Walker persistence exponent to second order 



As a first step in computing the mean and variance of S\ in the two ensembles we 
shall determine the structure of these quantities. For cither the restricted or the 
unrestricted average we have 

(5i)=S' f^' (31) 
4 id Jo 

in which the notation ^2' means the sum over combinations of the indices such 
that each index is different. (Note that all permutations of a set of distinct indices 
are included in this sum). Since we are averaging over ensembles of trajectories of 
independently diffusing particles, and the indices i and j are always different, we have 

^)=iE'r dr ' (i:j(T ' )>( ^ (r ' )> - (32) 

id Jo 

Hence to calculate the mean of Si we must find the mean velocity of a random walk 
with diffusion constant unity in each of the ensembles. 

For the variance of S\ we have 

(Si 2 ) - (5x> 2 = ^E'E' f dT i f dT s ({Mn)Mn)Mr2)Xe(r2)) 

id kj J ° Ja 

-{MTi)Mn))(X k fa)M'nt)j) ■ (33) 
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The terms in the double summation can be divided into three groups. In the first, all 
four indices are different, and two averages factorise to cancel. A second set of terms 
has two indices the same, the other two different (there are four ways this property 
can be satisfied). Finally there are two ways to arrange for the indices to comprise 
two pairs the same. Thus the variance can be written as 

((*>(n)Xi(T 2 » - (X,MHX,M)) (X i (n))(X i ( n )) 

i,j,k 

+2^' (/(XiMMTiMXMXjfo)) - (A > l (r 1 ))<X l (r 2 ))<A v J (r 1 ))(A v J (r 2 ))) (34) 

id 

which reveals that we must calculate not just the mean velocity {Xi(r')) but also the 
velocity correlation function (Xj(ri)X,(T2)) over the two ensembles. 

The unrestricted averages are easily found. Since each of the walkers performs diffusion 
with a diffusion constant of unity in the rescaled time r, we have Xi(r') = ^(t') where 
T)i(T f ) is a Gaussian white noise with zero mean and a correlator (??i(ri)?7i(r2)) = 
25(t2 — ti). Hence 

(X(t'))u = and {X 1 (t 1 )X 1 (t 2 )) u = 25{t 2 -t 1 ) . (35) 

The restricted averages involve more work, since we must calculate the mean velocity 
of a walker at time r' given that it survives until a later time r; similarly we must 
include the fact that the walker survives until a specified time r in the calculation 
of the velocity correlation function. Before getting into the details, therefore, let 
us propose their general form and the implications of such a form in context of the 
expression ([30|). On dimensional grounds, we suggest that 

(X i (r')) R = ^=F i (r /r , y/t) (36) 



(XiT^Xi^R = 2 



S(t2 - n) H Gi (t /t- ,t-/t+, t + /t) 

T_t_ 



(37) 



in which To is some early timescale in the problem and r_ (t + ) are the smaller 
(respectively larger) of t\ and r 2 . Guided by (|35|), and confirmed by explicit calculation 
to be presented below, we have included a delta-function contribution in the two-time 
correlation function. 



The upshot of this is that all the integrals in (|32|) and (|3_J) reduce to one of two forms. 
These are 

f Kf{ro/r'y/r) (38) 
Jo T 

f ^ f + dr_ g (T /T_,T_/T + , t+/t) (39) 
Jo T + Jo 

in which the functions / and g are various combinations of the velocity and velocity- 
correlation functions. It is not the values of these integrals that we are interested in, 
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but the limit that appears in the perturbation expansion (|30|). It is straightforward 
to show that 

lim r-f f ^/(tq/t', t'/t) = /(0, 0) (40) 
d /' r c!t f TJr 

Urn r— / / dr_ <?(t /t_, t_/t+, r+/r) =/ dw 5 (0,u,0) . (41) 

r^oo dr 7 t\ J Jo 

The fact that zeros appear in the arguments of the functions / and g implies that 
when calculating (Xi(r'))n and (Xj(ri)Xj*(T2)) r we need only determine their forms 
in the regime To <C t' , t\ , t% <C r to obtain the expansion ([30|) . 

We now indicate how to calculate these quantities. First, note that (X(t')) r = 
■£y(Xi(T')) r. This reduces the problem to the calculation of the mean position of a 
walker at time r' taking into account that it does not cross the origin (survives) until 
at least a time r. For the case where a walker's initial position Yi > 0, this quantity 
can be expressed as 

(X t (r'))n = fdX, X* PliT 'f;^ l) ( 42 ) 
Jo Qi{r\Yi) 

in which Pi(r; Xi,r'\Yi) is the probability that a single walker visits Xi at time r' and 
does not visit the origin at any time < r' < r. In order to average Xi with respect 
to this distribution, it has been normalised through the factor Qi{r\Yi), which is the 
probability the walker survives until time r and is given by equation (|li"|). 

Since the diffusion process is Markovian, we can perform the factorisation 
P x {T;Xi,T , \Y i ) = Qi(t - T , \X i )G 1 (X i ,T > \Yi), where we recall that G x {X U T , \Y i ) is 
that probability that a walker is at Xi at time r 1 given that it started at Yj, and 
whose form is given by equation (|^). Inserting the explicit expressions for Qi(t\Y) 
and Gipf,T|F) into © we find 



At' 



roo crf (ivfcO cxp ( - ~ exp 
(X i (r')) R = dXiXi L : li. (43) 



Now we can use the fact that when calculating the exponent 8 using ( |30| ) it is sufficient 
to know {Xi(r'))R in the regime To < r' < r. We take To = Yf (recalling that the 
rescaled time variable has units of length squared) and on expanding the integrand in 
the previous expression, we find 

W))r = I dX 4 X? exp (-^) = 4^- (44) 

in this regime. Finally, we differentiate to obtain 



(45) 
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when t <C t' <C t and > 0. Since we are working at times much larger than 
t = Y? , the dependence on the magnitude of Y t has dropped out of this expression. 
However, the sign of the average velocity clearly changes if we consider a walker with 
initial position Yi < 0. Hence the correct limiting form of the function F,(to/t',t'/t) 
in equation ( |36| ) is 

*i(0,0) = ^ (46) 
in which o~i — sgn(Y^). 

The two-point velocity correlation function (Xj(Ti)Xj(T2)).R is obtained in a similar 
way. We defer the details of the calculation to the appendix, quoting here only the 
result over the required range tq -C t\ , t-i <C t which is 



S{T 2 - Ti) + 



(47) 



Recall that t_ = min{Ti,T2} and t + = max{Ti,T2}. Note that this expression is in 
agreement with the proposal (|37| ) , and reveals the limiting form of the function G to 
be 



G i (0,r-/r + ,0) = -J^-l. (48) 



7T 



Once again, the dependence of this quantity on the walker's initial position Yi has 
dropped out; furthermore, the two-point correlation function is independent of the 
sign of Yi , so we can drop the subscript i on G. 

We now consider the first order term in (B(3). We have 



{Si) R -{S l ) u = - Y^ f dT ' (Mt'))r(Mt'))r 

= lE' / T ^(to/t',t'A)F,(t /t',t'A) . (49) 
id Jo T 

The unrestricted averages have disappeared from the right-hand side of this expression 
because (Xjjy = 0. Now, using equation (pf0|), we find that the limiting form of this 
expression required in the expansion ( ^p| ) reads 

lim tA ((S 1 ) r -(S 1 ) u ) =^'^(0,0)^(0,0) = ~Y, a i°i ■ ( 5 °) 

To perform the remaining summation, we note that aiUj is equal to +1 if walkers i,j 
both start on the same side of the origin, and is equal to — 1 if they start on opposite 
sides of the origin. If there are Nl (Nr) walkers initially to the left (right), there 
are Nl(Nl — 1) + Nr(Nr — 1) positive terms in the summation and 2N^Nr negative 
terms. Thus one finds that the coefficient of a^v in (|3^) is 

lim T |-«S 1 ) B -(5 1 > p ) = i(A a -JV) (51) 

T^oo dT 7T 
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Let us now turn to the second-order term in ( |30| ) . From (|^) we have the general form 
1 



((S 2 1 ) R -(S 1 ) 2 R -(S 2 1 ) U + (S 1 ) 2 U ) 



1 

32 



dn / dr ; 



4^ A ijk (T /T 1 ,T 1 ,T2,T 2 /T) + 2^ B V (WtI , Tl , T 2 , T 2 /t) 



i -J 



In the regime ro <Ti,T2 < t we have from (35), (Ma) and ( 

"7T(5(r2 -n) 2 r H 



Ay fc (0,T_,T + ,0) 

B tf (0,r_,T + ,0) = 



8(7 jCTfe 

16 



T2 



1-^-1 



(52) 

(53) 
(54) 



As with the first order term, we calculate the limit in ( |30| ) by making use of the results 
([fOj) and (|4l|). This procedure yields the result 

_d 1 

t^oo dr 2 



lim r^; «S?) R - {S 1 ) 2 R - {S!)u + {S 1 ) 2 U ) 



1 



7T + 4 / du- 



(55) 



To obtain the final expression for the expansion (pG) we use the fact that 



\\u ^~ 1 * 1 =2(ln2-l) 



(56) 



Also, we must consider the possible ways of choosing three particles labelled i,j,k 
from a set of N particles comprising a number Nl with a = — 1 and Nr with a = 1. 
These considerations lead us to 



J2'a j a k = (N-2)(A 2 -N) 



(57) 



Similarly the number of terms in the sum ■ is just N(N — 1) which finally allows 
us to write down an expression for the exponent 9 which constitutes the main result 
of this paper. It reads 

= E + I(iV - A 2 )a N 

i 7T 

+ 4j f(JV - 2)(JV - A 2 )(tt - 8(1 - ln2)) + 2N(N - 1)1 a 2 N . (58) 

7T 

As a check of this formula, we consider the exactly solvable cases of N = 1 and N = 2. 
When N = 1 , one can simply transform to a frame in which the particle is stationary. 
Then, its survival probability is given by equation ( |ll| ) with D replaced by D + D' 
which implies a persistence exponent of 9 = i. Since A 2 = 1 when N = 1, we find 
that the perturbative expression for 9 (158) also gives the exact result 9 = h. 
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When N = 2, there are two possible arrangements of the walkers. Either they surround 
the particle, in which case A = and the exact exponent is given by equation (|l7|), 
or they are both on one side of the particle, which has A = ±2 and 9 is given by 
equation (jlj). Expanding these exact expressions as series in a 2 one finds 

6 = 1 ± -a 2 + -U2 + • • • (59) 

where the plus sign is taken for the case A = and the minus sign for A = ±2. Thus 
to second order in 012 we find agreement between the perturbation series ( |58|) and the 
exactly known results. 



5. Discussion and conclusion 



In this work we have studied the survival probability of a diffusing particle in the 
presence of a finite number of mobile traps. In the absence of a complete exact 
solution for the problem, we have devised a method for calculating the persistence 
exponent 9 [defined by equation (|lj)] valid when the particle's diffusion constant is 
small. In order to calculate the coefficients in the expansion, it is necessary to find 
velocity correlation functions for a random walker in the presence of an absorbing 
boundary. In this paper, we calculated the exponent 9 to second order, culminating 
in the expression (p8|). We stress that, in principle, one could go to higher order, 
although the amount of work required to obtain n-time velocity correlation functions 
is likely to increase rapidly with n. 

We wish to make a few remarks about our result (pi). First, as we have already 



noted, the expansion parameter defined by equation (15) has an interpretation in 
the context of the (hyper)wedge geometry that has been previously used to obtain 
results for this problem. Although the connection appears at first sight intriguing, it 
is probably no coincidence since in the two approaches the parameter arises through a 
rescaling that gives rise to uncorrelated diffusion of the traps with diffusion constant 
unity. However we do see that a system with perturbed equations of motion (treated 
using the path-integral approach) can be related to one with the original equations of 
motion but perturbed boundaries (i.e. the hyperwedge geometry). 

Secondly, we find that the expansion to second order (|58| ) depends on two parameters: 
the total number of traps in the system N and the square difference A 2 in the number 
of traps initially placed on either side of the particle. The origin of this observation 
lies in the fact that at sufficiently large times the initial position of a particular trap 
is unimportant. Of course, the number of traps initially positioned to each side of the 
particle plays an important role, and enters into the perturbation expansion through 
summations in which some terms are positive and some negative depending on the 
initial arrangement of the traps. We would expect, therefore, the combination TV — A 2 
to appear in higher-order terms in the expansion. 



We also note from the expansion (58) that when the asymmetry is small, i.e. A 2 < N, 
the exponent 9 is an increasing function of the particle's diffusion constant D' at fixed 
trap diffusion constant D. Hence, in this near-symmetric situation, the particle is 
more likely to survive longer if it is at rest than if it is slowly diffusing. Of course, we 
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do not know if 9 continues to increase for large D' since we only have the first two 
terms in the perturbation expansion. However, our intuition suggests that this be the 
case. 

The case of large asymmetry, A 2 > N is more subtle since one of the terms in the 
coefficient of a 2 N is always positive. We learn a little by defining the critical asymmetry 
A c at fixed N and ajv to be that for which the persistence exponent 9 = y . Then for 
A > A c one has a regime in which the particle's survival probability is not (locally) 
maximised by remaining still. One can understand this phenomenon from the extreme 
case of all traps being on one side of the particle (A = ±JV). Then, there is some 
probability for the particle to diffuse into the region which is devoid of traps, thereby 
increasing its survival probability. By rearranging the expansion ( [58| ) one finds for 
small ajv that 



Of course, it is again difficult to make any statements for large a at given only the first 
two terms of the perturbation series. 

The final comment we wish to make about the form of ( |58| ) is that, although the 
quantity oln seems to be the natural expansion parameter in the problem, the presence 
of a nontrivial coefficient in the second order term is not suggestive of there being a 
simple, closed expression valid for all values of the parameters N, A, D and D' . 
Nevertheless, one might find that the marginal case N = A 2 is in some way much 
easier to understand. 

The trapping reaction studied in this work suggests a couple of generalisations. Firstly 
there is the question of how to treat higher dimensions. In a forthcoming work fl5f , 
we extend the methods introduced in JL4) for the case of an infinite sea of traps to 
dimensions greater than one. It would be interesting to see if the perturbation theory 
could be applied to higher dimensions to treat the case of a finite number of traps. 
Finally, we have studied in this work the trapping reaction on an infinite line. It 
appears that the problem in the periodic system is much harder to treat than the case 
considered here. 
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Appendix A. Two-point velocity function for surviving walks 

In this appendix we outline the derivation of the correlation function (X(ti)X(t2))r as 
given by equation (^) . Recall that this is the two-point velocity correlation function 
for a single random walker that has a diffusion constant Dq and does not cross the 
origin at any time < r' < r, where r > t%,T2- We shall assume, with no loss of 
generality, that the walker's initial position Y > and that t\ < t 2 . 




N + 



2N(N-1) 



a N + 0(a 2 N ) . 



(60) 
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The starting point is to consider the probability P\(r\ X\,T\ \ X 2 , t 2 \Y) for the walker 
to be at Xi at time t%, X 2 at time r 2 and to survive until time t given that it started 
at Y. Then, 



(X( n )X(T 2 )) R 



dX 1 / dX 2 X 1 X 2 



P 1 (t;X 1 ,t 1 ;X 2 ,t 2 \Y) 



dX x / dX 2 XiX 2 



Qx{t\Y) 

Qi{t - t 2 \X 2 )Gi(X 2 ,t 2 - nlx^dix^nlY) 



Qi(r\Y) 

in which the Markovian property of the diffusion process has been used to write the 
joint probability distribution function Pi(t; Xi, t\) X 2 , t 2 \Y) in terms of the functions 
G\(X,t\Y) and Qi(t\Y) defined through equations (^|) and (JTTj) . Note that it was 
also necessary to normalise the distribution according to the probability that a walker 
starting at Y survives a time r. 

Assuming that To = Y 2 -C t\,t 2 <C t, we can expand both the numerator and 
denominator in (AT) to first order in Y/y/ri and X 2 /^/t to find 

' ' dX 1 I dX 2 X 2 X^exp 



(X(n)X(T 2 )) R 



4:1TTly/Ti(T 2 - T\) JO 
(X 2 - Xl) 2 

4(r 2 - ti) 



Xf 
4rj 



exp 



exp 



(X 2 + X 1 ) 2 
4(r 2 -ri) 



(A.2) 



This integral can be evaluated without recourse to further approximations. First one 
integrates over X 2 , which yields 



(X( Ti )X(t 2 )) 



Ft 



x erfc 



V4(r 2 -n) 



4(T2 - Ti) 



(^ + 2(r 2 - n )) 



Ai exp I — 



4(r 2 - n) 



(A.3) 



To perform the remaining integration, it is helpful to make a change of variable 
Xi = 2^/t 2 — T\u which leads to 



(X(n)X(r 2 )) R = 



16(r 2 -n) 



5/2 fco 



du 



\/7T 2 



(A.4) 



in which a = t 2 /t\ and f3 = (r 2 — t\)/t\. Then, using 
1 

~2rf 



duu 3 e- au 

duu^e-^eviu 
we obtain, after a lot of manipulation 

4ri 



d 

"dp 



due ^ u eriu 



d_ 

"d/3 



arctan(l/v^) 



(X(n)X(r 2 )) 



R 



Tl 

n 



arctan 



t 2 - Tl 



T 2 - Tl 



(A.5) 



(A.6) 



(A.7) 
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To find the velocity correlation function (X(ti)X(t 2 ))r we differentiate with respect 
to Ti and t 2 . We find 



4 f arctan./Tv + ^F^ n < r 2 

dT2 n { 2 arctan + 2^-^ n > t 2 

For Ti < t 2 one has 



(X(n)X(r 2 )) R = — J r ± (A.9) 

7TT 2 V Tl 

and the same expression with ti <-> t 2 when ti > t 2 . At ti = t 2 , however, there is a 



jump in (A. 8) of height 2 which implies that 

(X( n )X(T 2 ))R = 2 



2 

<H7"2 - Tl) H 



7TT_|_ y T_ 



(A.IO) 



Thus we conclude our derivation of (47) 
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